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^ ; ABSTRACT 

■ Context. The electron-cyclotron maser instability is widespread in the Universe, producing, e.g., radio emission of the magnetized 

rS^ ■ planets and cool substellar objects. Diagnosing the parameters of astrophysical radio sources requires comprehensive nonlinear sim- 

ulations of the radiation process taking into account the source geometry, 
ly-^ \ Aims. We simulate the electron-cyclotron maser instability (i.e., the amplification of electromagnetic waves and the relaxation of an 

unstable electron distribution) in a very low-beta plasma. The model used takes into account the radiation escape from the source 
region and the particle flow through this region. 

Methods. We developed a kinetic code to simulate the time evolution of an electron distribution in a radio emission source. The 
\ model includes the terms describing the particle injection to and escape from the emission source region. The spatial escape of the 

. , emission from the source is taken into account by using a finite amplification time. The unstable electron distribution of the horseshoe 

• type is considered. A number of simulations were performed for different parameter sets typical of the magnetospheres of planets and 

{Zlj ultracool dwarfs. 

Results. The generated emission (corresponding to the fundamental extraordinary mode) has a frequency close to the electron cy- 
clotron frequency and propagates across the magnetic field. Shortly after the onset of a simulation, the electron distribution reaches 
a quasi-stationary state. If the emission source region is relatively small, the resulting electron distribution is similar to that of the 



o 



^ ■ injected electrons and the emission intensity is low. In larger sources, the electron distribution may become nearly flat due to the 

' wave-particle interaction, while the conversion efficiency of the particle energy flux into waves reaches 10 - 20%. We found good 

agreement of our model with the in situ observations in the source regions of auroral radio emissions of the Earth and Saturn. The 
expected characteristics of the electron distributions in the magnetospheres of ultracool dwarfs were obtained. 
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1 . Introduction The electron-cyclotron maser emission is highly sensitive to 

the parameters of plasma, magnetic field, and energetic particles 



O 

■ in its source. However, utilizing the high diagnostic potential of 
^ ■ The electron-cyclotron maser instability (ECMI) is believed to this emission is a difficult task. Although some parameters of the 

be responsible for generation of auroral radi o emissions from emission source can be estimated using general features of the 

. ^ . the magnetized planets of th e Solar System (|Wu & Lee| |l979t maser mechanism (e.g., the emission frequency, as a rule, is very 

; |Zarka||1998|; |Treumann' 2006). It is very likely that similar ra- close to the electi-on cyclotron frequency which allows us to find 

■ dio emissions are also typical o f exoplanets ( Bastian et al. ,2000; the magnetic field strength), most of the source parameters affect 
C3 ; |ZarkaetalJ|200H; |Zarka| |2007|; |Hess & Zarka||201l|), although the emission properties in a compHcated and nonlinear way. A 



the sensitivity of the existing instruments does not allow us to de- comprehensive nonlinear simulation of the ECMI requires con- 

tect such radio signals. It was recently discovered that a number sidering the spatial movement of the waves and particles, which 

of very low-mass stars and brown dwarfs (collectively known as requires a lot of computational resources, 
ultracool dwarfs) are the sour ces of intense periodi c radio bursts 

in the GHz frequency ra nge (iHallinan et alJl2006l l2007l l2008t In most works devoted to nonlinear simulation of the 

iBerger et al.ll2008ll2009 h: these bursts seem to be produced due ECMI, the spatial movement of waves and particles was sim- 

to the ECMI in a way similar to the planetary auror al radio e mis- ply ignored, which corresponds to an infinite and homoge- 

sions , but in a much stronger magnetic field (IKuznetsov et all neous source (the so-called diffusive limit). Kinetic simulations 

l2012h . The electron-cyclotron maser has been also suggested of relaxation of the loss-cone e lectr on distribution in a dif - 



as the generation mechanism of the periodic radio bursts frorn fusive hmit were performed by [Asch wanden & Ben3 (Il988l) : 

the magnetic Ap star CU Virginis (TrigiHo et al.' '2000l. 1201 it lAschwandenl(ll99ot) : rFIeishman & Arzneii(i2000l) . Similar simu- 

|Ravi et al. 2010) and of cert ain types of radio bursts occurring lations f or the horse s hoe-li ke distribution are presented in the 

during solar and stellar flares dMelrose & Dulkl 1 982l : iDuIkl 1 9851: paper of iKuznetsovl (1201 Ih (hereafter Paper I). A number of 

[Fleishman & MeInikovlll998h . simulations in a similar approach (i.e., ignoring the source 
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inhomogeneity and finite size) for different electron distribu- 
tions were performed using the particle-in-cell technique (e.g., 
[Pritchett 1984a b, 1986b; Pritchett & Strangeway 1985). More 
complicated particle-in-cell simulations considered an inhomo- 
geneous emission source and used specific boundary conditions 
and/or numerical methods to account for the radiation escape 
from the source as well as for the particle flow through the 
emission source region (.Pritchett 1986a : Wi nglee et al.. ,1988; 
iPritchett & Wingleelll989tlPritchett et aLli2002l) . 



In this work, we use a different approach. Like in Paper I, the 
particles and waves are simulated using a kinetic description. 
The used model does not explicitly include the spatial move- 
ment of the particles and waves and their spatial ditributions 
within the emission source, i.e., the parameters of the particles 
and waves are treated as averaged over the source volume. The 
spatial movement of the electromagnetic waves (which results 
in their escape from the source) is taken into account implic- 
itly by introducing a finite amplification time. In turn, the spatial 
movement of the electrons is taken into account implicitly by 
introducing the injection and escape terms in the kinetic equa- 
tion. Although the considered model is still oversimplified, we 
believe it is a much better approximation to reality than the 
diffusive limit. Like in Paper I, we consider the horseshoe-like 
distribution of the energetic electrons (e.g., Delorv et al. 1998; 
Ifergun et al.ll2000t : iStrangewav et aLll2001l) . The dispersion rela- 
tion of the electromagnetic waves is assumed to be similar to that 
in a vacuum. Simulations are performed for different parame- 
ter sets corresponding to the plasma parameters measured in the 
sources of the auroral radio emissions of the Earth and Saturn, 
and to the expected parameters of the magnetospheres of ultra- 
cool dwarfs. Although our model, in principle, can be applied 
to non-stationary processes, we focus in this work on finding 
quasi-stationary solutions and analyzing their characteristics. 

Generation of terrestrial auroral kilometric radiati on in fi- 
nite-si z e low - density plasma cavities was stu died by ICalverll 
(1981', '1982); Louarn & Le Oueau| (Il996allbl) : iLouarn' (2006*); 
[Burinskava & Rauch (2007), who focused on the wave disper- 
sion and propagation and the conditions of radiation escape from 
the cavity, but did not consider the nonlinear wave-particle in- 
teraction. The difference of this work is that, firstly, our aim is 
to investigate the nonlinear wave-particle interaction, while the 
wave propagation effects are assumed to be incorporated into a 
single factor - the wave amplification time. Secondly, we do not 
consider the possible formation of a discrete spectrum of eigen- 
modes in a plasma cavity, assuming that the quality factor of 
such a resonator is too low due to strong energy leakage. 

As can be seen from above, we do not consider formation 
of fine temporal and spectral structures in the dynamic spec- 
tra of radio emission. It was recently proposed that these fea- 
tures may be a reflection of small-scale short-lived nonlinear 
structures (such as electron holes) in the emission source (e.g., 
iTreumann et akllToi ll |2012|) . However, simulating the electron- 
cyclotron maser emission from (or in presence of) electron holes 
requires a more sophisticated approach and thus is beyond the 
scope of this paper 

The numerical model used in the simulations is largely simi- 
lar to that in Paper I; therefore, most formulae and technical de- 
tails are omitted and only some important facts and expressions 
are briefly repeated. The used model is described in Section |2] 
The simulation results (for a wide range of parameters) are pre- 
sented in Section |3] The comparison of the simulation results 
with the observations is made in Section |4] The conclusions are 
drawn in Section|5] 



2. Model 

2.1. Particle dynamics 

In the general case, the time evolution of the electron distribu- 
tion in a stellar/planetary magnetosphere should be described 
by the Fokker-Planck equation including the effects of the mag- 
netic field convergence, magnetic-field-aligned electric field, and 
wave-particle interaction (e.g.. iHamilton et al..il990) . However, 
as said above, solving the exact equation would require enor- 
mous computational resources due to necessity to consider ex- 
plicitly the spatial movement of the particles. In this work, we 
use a simplified approach where the spatial dependence of the 
electron distribution function is neglected and the particle move- 
ment is treated approximately using the particle fluxes in and out 
the considered volume. Therefore, the kinetic equation for the 
distribution function / can be written in the form 



dt 



5/(p,f) 



dt 



5/(p,f) 



dt 



dfipj) 



dt 



(1) 



where p is the electron momentum and the terms (5//(9f)inj, 
(5//(9f)esc, and (c?//i9f)rei describe the injection of the energetic 
electrons into the considered volume, escape of the electrons 
from this volume, and variation of the distribution function due 
to the wave-particle interactions, respectively. The effects of the 
converging magnetic field and the parallel electric field are as- 
sumed to be incorporated into the terms (df/ (9r)inj and {df/dt)^^^, 
since the escaped electrons can be then reflected from the mag- 
netic mirror and/or re-accelerated by the electric field, so that 
they re-appear in the simulation model as the injected electrons 
(but with a different value of the momentum p). 

We assume that the particle injection is time-independent, so 
that the corresponding term in Eq. ([T) can be written in the form 



\dmt)] 
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dt 


inj 


\dt A 



(2) 
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where [drif^/dt^nj is the constant injection rate and is the nor- 
malized to unity distribution function of the injected particles. 

The particle escape from the radio emission source can be 
described by the expression 



dfipj) 



dt 



/(P,0 

Tesc(p)' 



(3) 



where Tesc is the characteristic escape time. For the particles with 
constant (time-independent) velocities, the escape time can be 
estimated as Tesc - Rz/l^zl where is the longitudinal (i.e., 
along the magnetic field) source size and u~ is the longitudi- 
nal component of the velocity. However, in situ measurements 
within the sources of terrestrial and Saturnian auroral radio 
emissions reveal the presence of magn etic-field-aligned electric 
fields which accel erate the electr ons (Ergun et al." 19981, 120001; 
Treumann 2006; L amv et aP 1201 0; Schippers et al. 201 1]); as a 
result, V, + const and even the particles with zero initial velocity 
eventually leave the source region. For simplicity, we assume in 
this work that the particle escape time is constant and is given by 
the expression 



Tesc ^ R:/Vb, 



(4) 



where is a typical particle speed. In a quasi-stationary state 
(when the electron density is constant), the described model is 
equivalent to t he "recycling" metho d used in the particle-in-cell 
simulations bv lWinglee et al.l (Il988h . 
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2.2. Wave dynamics 

If spatial movement of waves is neglected, the kinetic equation 
for the waves of mode cr can be written in the form 



It 



/'^'(k,f)w<">(k,o, 



(5) 



where W^*^' is the energy density of the waves in the space of 
wave vectors, k is the wave vector, and is the growth rate 
(i.e., the above equation describes wave amplification with time 
due to an instability). The wave-mode index cr is hereafter omit- 
ted for brevity. 

In Paper I, we considered a diffusive limit, when the wave 
amplification was limited only by nonlinear effects - relax- 
ation of an unstable electron distribution due to interaction with 
waves, which resulted in a decrease of the growth rate. In this 
work, we consider an opposite limit, when the waves are not 
accumulated in the emission source. Escape of the waves from 
the amplification region (i.e., from the volume occupied by the 
electrons with an unstable distribution) naturally limits their in- 
tensity. To account for this effect, we assume that (i) the elec- 
tromagnetic waves are amplified due to the electron-cyclotron 
maser instability only during a small time interval Af and (ii) the 
growth rate remains nearly constant during this time interval. As 
a result, the maximum energy density of the waves can be esti- 
mated as 



^^'"(k,?) 



Wo(k)exp[r(k,f)Af(k,f)]: 



(6) 



where Wo is some initial energy density (which is determined, 
e.g., by the spontaneous emission processes). 

The wave amplification time can be estimated as 



Af ^ min 



drjdt 



dr^/dt 



dkjdt 



dk^/dt 



(7) 



where R- _l are the characteristic sizes of the emission source 
(along and across the magnetic field, respectively), Ak,^j_ are the 
characteristic sizes of the wave amplification region in the space 
of wave vectors, and dr^^^ldt and dk^ j^/dt are the corresponding 
movement speeds of a wave packet. Escape of the waves from 
resonance with energetic electrons due to variation of their wave 
vector in an inhomogeneous medium can be an important fac- 
tor affecting, e.g., formation of some types of solar radio bursts 
(IVlasov et alJl2002l: iKuznetsovlIlM l2008h . However, we con- 
sider in this work the case when the plasma density is very low 
so that the "vacuum approximation" for the wave dispersion can 
be used (see below) and the wave vector is nearly independent 
on the medium parameters and thus on the spatial coordinates; 
as a result, 5A: ,x — > and the wave escape from resonance in 
the space of wave vectors is negligible. Also, we consider ring- 
like or horseshoe-like electron distributions that generate waves 
mainly in the direction perpendicular to the magnetic field (see, 
e.g.. Paper I). As a result, \dr,/dt\ <s: \drj^/dt\ and only the spatial 
escape of the emission from the source in the transverse direc- 
tion needs to be considered. The expression for the amphfication 
time is then reduced to 



Af 



sin 6 



(8) 



where Vg^ is the group speed of the waves and is their propa- 
gation angle relative to the magnetic field. We should note that, 
according to the satellite measurements, the terrestrial auroral 
kilometric radiation is generated in auroral cavities surrounded 



by a denser plasma (e.g.. lErgun et al.l[T998Ll200(j|) . The electro- 
magnetic waves can be reflected from the walls of these cavities 
and return back thus increasing the time spent in the emission 
source. However, the reflected waves will have a different (not 
perpendicular to the magnetic field) direction of propagation and 
thus will not further participate in resonance with the electrons. 
Therefore the wave amplification time will still be determined 
by Eq. ([8]), which now will characterize the time of wave escape 
from the amplification region in the space of wave vectors. 

The formula ^ represents the energy density of the waves 
leaving the amplification region and thus determines the emis- 
sion intensity outside the source. By using Eq. (|5]), we can find 
the total energy radiated from unit volume per unit time: 



dW(t) 



dt 



rad 



= J r(k,f)w^^''^(k,f)d3 



(9) 



2.3. Wave-particle interaction 

Variation of the electron distribution due to interaction with 
the electromagneti c waves is described by the equation (e.g., 
lAschwandenll 19901) 



\df(p,t)] 
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dt 


rel 


dpi 



Dijipj) 



5/(p,f) 



dpj 



where Dij is the diffusion tensor: 
A/p,f) = J 



^^^'^\ikjW^(k,t)d\, 



(10) 



(11) 



o) is the wave frequency and ^ is the factor describing the effi- 
ciency of the wave-particle interaction (it includes the resonance 
condition and summation over cyclotron harmonics). If several 
wave modes can exist simultaneously, then the diffusion tensor 
Dij is the sum of the tensors ( fTTT i for the different modes. In turn, 
the growth rate of the electromagnetic waves due to intera ction 
with the energetic electrons equals (e.g.. lAschwandenll 1 990l) 



-I 



r(k,f)= Ji{Kp)ki 



^dV 

dpi 



(12) 



The exact expressions for the factor and the methods of cal- 
culating the growth rate, diffusion tenso r coefficients, and diffu- 
sion rate are g iven in many works (e.g.. 'Melro se & Dulklll982l : 
Asc hwanden & Benz 1988; Aschwanden 1990) and summarized 
in Paper I. 

When calculating the diffusion tensor (fTTT i. we assume that 
the wave energy density Wk equals its maximum value W^™"' 
given by Eq. since the waves with the highest energy density 
make the maximum contribution into the relaxation process of 
an unstable electron distribution. Thus, by using all the assump- 
tions listed in Sections l2TT]|2.3l the time derivative of the electron 
distribution function df/dt is reduced to an explicit function of 
the distribution function itself; in turn, the simulation model is 
reduced to the differential equation ([Tj). In the below simulations, 
the electron distribution function /(p) is defined on a regular grid 
in (p, a) space, where a is the electron pitch angle. Equation ^ 
is numerically integrated with respect to time using the Gear for- 
mulae of fourth order with an adaptive stepsize; the details of the 
numerical code are given in Paper I. 

2.4. Model parameters 

In situ measurements within the sources of terrestrial and Sa- 
turnian auroral radio emissions have revealed that such regions 
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have the followin g characteristic features (e.g.. lTreumanDl2006t 
iLamv et a l?'2010'): (i) the plasma density is very low so that the 
electron plasma frequency ojp is much lower than the electron 
cyclotron frequency wb ; (ii) the cold (thermal) electron compo- 
nent is almost absent, so that the energetic electrons with the 
energy of a few keV dominate. Under these conditions, the dis- 
persion of the electromagnetic waves differs significantly from 
that in a cold plasma (Pritchett 1984a b; Le Oueau et al. 1984a,b; 
Winglee 1985; Strangeway 1985, 1986; Robinson 1986 J 19871: 
Le O ueau & Louam .1989: Louarn & Le Oueau 1996b). In par- 
ticular, numerical simulations by Robinson (1986, 1987) have 
shown that if the plasma parameters satisfy the condition 



(13) 



then the dispersion of the waves (in the case of quasi-transverse 
propagation) becomes similar to that in a vacuum, i.e., the re- 
fraction index becomes close to unity for both the ordinary 
and extraordinary waves, the low-frequency cutoff of the X- 
mode (which is predicted by the cold plasma theory) disap- 
pears and the dispersion branches corresponding to the X- and 
Z-modes of a cold plasma reconnect forming a single branch 
(with only a small wiggle near the electron cyclotron frequency). 
Condition ( fT3l l is well satisfied within the sourc es of terres- 
trial and Saturnian auro r al radio emission s (e.g., iDelorv et alj 
Il998t lErgun etalJlIOOOl: iLamv et alJlIoTol) : most likely, radio 
emissions of other magnetized planets and ultracool dwarfs are 
produced under similar conditions. In addition, observations 
of the waves within the sources of terrestrial auroral emission 
have found n o indication of the low-frequency X-mode cutoff 
dErgun et al.i ri998). Therefore, we adopt in this work a vacuum- 
like dispersion relation, with N = 1 and Vg^ - c for both the 
ordinary and extraordinary modes. Other details of calculating 
the growth rate and diffusion tensor for such waves are given in 
Paper I. 

We assume in this work that the inje cted energetic elec- 
trons have a horsesh oe-like di stribution dPelorv et al.l Il998t 
lErgun et ani2000: Strangewav et a l. 2001). which is modeled by 
the following expression (the same as in Paper I): 



/inj(p) = Aexp 



(p-Phf 



Ap 




(14) 



where p\, and Apb are the typical momentum of the injected elec- 
trons and their dispersion in momentum, respectively, fi - cos a, 
fic - cos ffc, Q-c and Ap^ are the loss-cone boundary and the 
boundary width, respectively, and A is the normalization factor 
The electrons with < a < 7r/2 move upwards (i.e., toward a 
decreasing magnetic field), and the electrons with 7r/2 < a < n 
move downwards (toward an increasing magnetic field). If - 
0, we obtain an isotropic ring-like distribution. 

The initial wave energy density Wo depends on several fac- 
tors, including the spontaneous radiation and the emission com- 
ing from beyond the considered region. We characterize this en- 
ergy density by the corresponding effective temperature Tq, i.e. 



Wo 



(27r)3 



const. 



(15) 



where is the Boltzmann constant. 

When solving the differential equation ([TJ, we are interested 
in finding a quasi-stationary state with df/dt — » 0. Therefore, 
the initial electron distribution function (at f = 0) is unimportant. 
We either assumed that the initial distribution function was equal 



to zero, or used a result of another simulation as the initial condi- 
tion for the next simulation run, in order to provide faster conver- 
gence. Simulations were stopped when the variation rate of the 
distribution function fell below a certain threshold (10^^ - 10"^ 
of the injection rate). 



3. Simulation results 

3.1. Electron distributions and emission cliaracteristics in a 
quasi-stationary state 

It is easy to show analytically (and is confirmed by our numerical 
simulations) that for the above described model, with an increas- 
ing time, the electron distribution asymptotically approaches a 
quasi-stationary state with df/dt —> regardless of the initial 
distribution, with Tesc being the characteristic timescale of this 
process. The same conclusion can be made for the emission in- 
tensity and spectrum, since they are determined by the electron 
distribution at a given time. Thus if the characteristic variation 
timescales of the observed radio emission far exceed the particle 
escape time Tesc (which implies that the source parameters are 
stable at those timescales), we can expect that the electron dis- 
tribution in the emission source has reached a quasi-stationary 
state or is very close to it. Radio emissions of the planets and 
ultracool dwarfs, as a rule, satisfy this condition (see Section|4]i. 
Therefore we focus in this work on the electron distributions and 
the emission intensities and spectra in a quasi-stationary state. 

Firstly, we consider the relative effect of the two factors af- 
fecting the electron distribution: particle escape from the radio 
emission source and relaxation of the electron distribution due to 
wave-particle interaction. We use the following simulation pa- 
rameters (that were chosen so as to reproduce the observed char- 
acteristics of radio emission from ultracool dwarfs, see Section 
IHfor details): the electron cyclotron frequency vb = 4.5 GHz, 
the transverse size of the emission source - 1000 km, the 
initial effective temperature of the waves To - 10^ K, the typical 
energy of the injected electrons E\, - 10 keV, the relative disper- 
sion of the electrons in momentum Apb/ pb = 0.2, the loss-cone 
boundary - 60°, the loss-cone boundary width AyUc = 0.2, 
and the electron injection rate (c?ne/i9f)mj = 5 x 10* cm""^ s"'. 
The particle escape time Tesc is variable. Since the extraordi- 
nary mode has been found to be strongly dominant under the 
above conditions (see Section [J!2] |. we consider only this mode 
in most of our simulations. Figure [1] shows the quasi-stationary 
electron distributions corresponding to the different values of 
the particle escape time. Note that the distribution function in 
each panel of Fig. [T]is normalized to its maximum value, while 
the actual values of the distribution function increase with Tesc 
which is indicated by the increasing total electron density n^. 
Figure |2] shows the coiTesponding growth rates of the extraordi- 
nary waves. Finally, Figure |3] summarizes the simulation results 
for the different escape times and shows the maximum growth 
rate and the total radiation energy flux from unit volume calcu- 
lated using Eq. (|9]i- The parameter A in Figs. |2][3] is the wave 
amplification factor that is defined by A = Wk/Wo; according to 
Eq. ©, InA = yAt. 

To avoid a possible confusion, we highlight here that Fig. 
[T]does not represent a temporal evolution of an electron distri- 
bution. Instead, each panel in Fig. [T] represents a final (quasi- 
stationary) stage of such evolution. These final distributions 
coiTespond to different operating conditions of the electron- 
cyclotron maser, namely, to different longitudinal source sizes R~ 
which results in different values of the particle escape time Tesc- 
Accordingly, Fig.|2]represents the final (quasi-stationary) growth 
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Fig. 1. Quasi-stationary electron distributions for the different values of the particle escape time. Each panel represents a normalized 
distribution function //max(/) vs. the normalized momentum pl{mgc) and the pitch-angle a. The simulation parameters are given 
in Section im 



rates, while Figs.[3]|5]plot some characteristic parameters of the 
quasi-stationary solutions vs. the emission source parameter resc- 
Figures [T^ and correspond to a case when the timescale 
Tesc is very small and thus the particle escape from the emis- 
sion source region is very fast. In this case, the electron den- 
sity (and, consequently, the wave growth rate) cannot reach a 
level sufficient to provide a significant wave amplification. As a 
result, relaxation of the electron distribution due to interaction 
with the waves is negligible and the electron distribution in a 
quasi-stationary state does not differ from the distribution of the 
injected electrons . The corresponding electron density equals 



i dt 



(16) 



The growth rate has a sharp peak at ^ 90° and oj/ais - 0.985 
that corresponds to the relativistic cyclotron frequency of the 
electrons with the energy of about 10 keV. The emission inten- 
sity is very low (see Fig. [3]). 

The quasi-stationary growth rate of the waves increases lin- 
early with Tesc, until the corresponding wave energy density ex- 
ceeds a certain threshold so that the wave-particle interaction 
begins to affect the electron distribution. Under the considered 
conditions, this happens at In A > 19. Figure[TJ) shows an exam- 
ple of a weakly-relaxed quasi-stationary state when the electron 
distribution is only slightly distorted in comparison with the dis- 
tribution of the injected electrons. In turn, the relaxation of the 



electron distribution reduces the wave growth rate, so that the 
growth rate peak broadens and flattens (see Fig. |2}5). Note that 
the wave energy density depends exponentially on the growth 
rate, and therefore only the waves in a very narrow spectral 
(A(l>/cl> ^ 0.01) and angular (A9 - 3°) range actually make a 
significant contribution into the total emission intensity as well 
as into the relaxation of the electron distribution. Although the 
fraction of the particle energy flux that is transferred to waves 
is relatively low (about 1% for the state shown in Figs. [T}? and 
I2J3), this conversion efficiency can be sufficient to provide, e.g., 
the observed intensity of terrestrial and Saturnian auroral kilo- 
metric radio emissions ( Gurnett 1974; Benson & Calvert 197^ 
iKurth et alJl2005HLamv et alJl201 Ih . 

For even larger values of the timescale Tesc (Figs. [T];-[T]J), 
the electrons interact with the waves for a longer time and thus 
drift further in velocity space (towards lower velocities) before 
they eventually leave the radio emission source. As a result, the 
particle dispersion in velocity increases with Tesc; for the case 
shown in Fig. [TJi, the empty space inside the "horseshoe" is al- 
most filled. The maximum growth rate remains nearly constant, 
but the flattened region around the peak of the growth rate broad- 
ens (Figs. 12}; -IJ}!); therefore, the total emission intensity gradu- 
ally increases with Tesc- 

Figures [1^ and correspond to the value of the timescale 
Tesc that provides the maximum emission intensity; the conver- 
sion efficiency of the particle energy flux into waves in the quasi- 
stationary state reaches 14.5%. One can see that the empty space 
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Fig. 2. Quasi-stationary growth rates of the extraordinary waves for the different values of the particle escape time. Each panel 
represents a growth rate y vs. the normalized frequency w/wb and the propagation direction 6. The simulation parameters are given 
in Section |37T] 




"'"esc S 



Fig. 3. Maximum growth rate (top) and total emission intensity 
(bottom) of the extraordinary waves in a quasi-stationary state 
vs. the particle escape time. The vertical lines correspond to the 
values of Tesc shown in Figs.[T]|2] The simulation parameters are 
given in Section lTTI 



around t; = is now completely filled by the electrons; the 
quasi-stationary electron distribution does not look like a "horse- 
shoe" but rather represents a combination of a nearly flat plateau 
and a downgoing electron beam (with energy of ~ 10 keV and 




10"^ 10"^ 10"^ 10"^ 10"' 



Fig. 4. Maximum growth rate of different wave modes in a quasi- 
stationary state vs. the particle escape time. The simulation pa- 
rameters are given in Section lTTl 

the pitch-angle dispersion of ~ 45°). With a further increase of 
Tesc, the electron distribution does not change qualitatively - the 
plateau level steadily rises while the beam-like component re- 
mains nearly the same (see Fig.[TJ). The peak of the wave growth 
rate becomes narrower (Fig. |2]F) and the radio emission char- 
acteristics become similar to those for a beam-driven electron- 
cyclotron instability, i.e., the emission becomes directed slightly 
downwards (with the maximum at ^ 90° - 93°). Although 
the maximum growth rate slightly increases with Tesc, the to- 
tal emission intensity steadily decreases but remains sufficiently 
high (with an energy conversion efficiency of > 10%). 

Equation ( fTSI l for the electron density in a quasi-stationary 
state Woo was obtained under the assumption that the wave- 
particle interaction is negligible. If such interaction is important. 
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Fig. 5. Maximum growth rates and total emission intensities of the extraordinary waves in a quasi-stationary state vs. the particle 
escape time for the different parameter sets. In each panel, most simulation parameters are the same as in Section 13.11 and one 
parameter (whose values are specified in the panel) is variable. 



the quasi-stationary electron density slightly exceeds Hoc (but by 
no more than 10"^). 

Note that the values of the particle escape time used in 
this Section were chosen rather arbitrarily. We have found that 
different ratios of the particle injection and escape rates, to- 
gether with the wave-particle interaction, yield qualitatively dif- 
ferent electron distributions and radio emission characteristics. 
Consequently, the actual characteristics of the particles and ra- 
dio emission in a planetary or stellar magnetosphere will be de- 
termined by the local conditions (such as the emission source 
dimensions, electron acceleration rate, etc.). In particular, the 
electron distributions observed in terrestrial or Saturnian mag- 
netosphere seem to be weakly or moderately relaxed like those 
shown in Figs.[TJ)-[T];, while the electron distribution in the mag- 
netospheres of ultracool dwarfs are expected to be strongly re- 
laxed like that shown in Fig. [If (see Section|4]for details). 



Figure |4] demonstrates the quasi-stationary growth rates of 
the considered wave modes corresponding to the different val- 
ues of the particle escape time Tesc- One can see that the ratios 
of the growth rates are nearly independent on Tesc, i e-, they are 
nearly the same for both the weakly and strongly relaxed elec- 
tron distributions. Although the maximum growth rate of the X2- 
mode somewhat increases at large Tesc, it always remains much 
lower (by more than an order of magnitude) than that of the XI- 
mode. Since the emission intensity depends exponentially on the 
growth rate, the X2-mode intensity is lower than that of the XI- 
mode by a factor of ~ 10'^. The growth rates and hence intensities 
of the Ol- and 02-modes are even lower. Thus we conclude that 
the XI -mode is strongly dominant under the considered condi- 
tions. The intensities of other modes are negligible, and therefore 
they do not make a contribution into both the escaping radiation 
and the relaxation of the electron distribution. 



3.2. Relative contribution of different wave modes 



3.3. Effects of varying ttie model parameters 



In a low-density plasma (with Wp/we ^ 1), the funda- 
mental extraordinary mode is expected to be the dominant 



mode of the ele ctron-cycl otron maser (e.g., iWu & Lei 



Fleishman & Yastrebovi 1 19941 : iKuznetsovl 201 ih . This conclu 



Melrose & Duii3ri9 82: Sha rma & VlahosiT 984t IWinglei 



1979; 



1985 



sion is generally supported by the o bserva t ions of the plan- 
etary auroral rad i o emissions dZarkal 11998 1 : iTreumannI l2006t 
lErgun et al J [20001 iLamv et alJl201 ij). In order to investigate a 
possible contribution of other wave modes, we performed a 
simulation considering four modes simultaneously: fundamental 
and harmonic extraordinary modes (XI and X2), and fundamen- 
tal and harmonic ordinary modes (Ol and 02). The simulation 
parameters were the same as in the previous Section. 



In Section 13.11 we considered the particle and emission char- 
acteristics for the different values of the particle escape time 
from the emission source, i.e., for the different longitudinal sizes 
of the source. Now we consider the effect of varying the other 
model parameters. Figure [5] is a collection of plots similar to 
those in Fig. [3] but corresponding to the different parameter sets. 
Each panel of Fig.[5]demonstrates the effect of varying a partic- 
ular model parameter, while the other parameters are assumed to 
be the same as in Section ITTI 

Figure [5^ demonstrates the effect of varying the particle in- 
jection rate {dn^ldt)my As expected, an increase of this param- 
eter results in an increase of the emission intensity; the quasi- 
stationary growth rate of the waves also somewhat increases. In 
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addition, a higher particle injection rate allows the emission to 
be efficiently produced in smaller sources (i.e., with a lower par- 
ticle escape timescale Tesc)- 

Figure |5j> demonstrates the effect of varying the transverse 
size of the emission source R± and hence the wave amplifica- 
tion time Af. If the particle escape timescale is relatively low, an 
increase of the transverse source size makes the wave amplifica- 
tion more efficient and thus increases the emission intensity sig- 
nificantly. In contrast, if Tesc is large, a larger transverse source 
size means a stronger relaxation of the electron distribution; this 
may result in a slight decrease of the emission intensity. Note 
that we talk here about the average emission energy flux from 
unit volume (5W/5f)iad; the emission intensity from the entire 
source volume is proportional to (dW/dt\adR]_ and therefore al- 
ways increases with Rj^. Also note that the actual parameter af- 
fecting the wave energy density is the amplification time Af. If 
the wave dispersion differs from the "vacuum approximation" 
used in this work, the amplification time may be also different. 
However, the effect of increasing or decreasing the wave ampli- 
fication time will be the same as that shown in Fig.jSj^. 

Figure|5}; demonstrates the effect of varying the initial effec- 
tive temperature of the waves Tq. One can see that this parame- 
ter affects the quasi-stationary growth rate of the waves provided 
that the wave energy density has reached the relaxation thresh- 
old (i.e., it is sufficiently high to affect the electron distribution). 
A higher initial intensity means that the waves need less ampli- 
fication to reach the mentioned threshold. On the other hand, the 
emission intensity itself is almost independent on Tq. Therefore 
we conclude that the initial effective temperature of the waves 
may be chosen rather arbitrarily; in most simulations, we set it 
to 10^ K. 

Figure |5]l shows the simulation results for the different en- 
ergies of the injected electrons E\,- Since the particle injection 
rate (5ne/5f)inj is assumed to be constant, an increase of the en- 
ergy Eb corresponds to an increase of the energy injection rate. 
In turn, this results in an increase of the emission intensity. The 
maximum conversion efficiency of the particle energy flux into 
waves is about 14.5% for all considered values of fib- 
Figure |5^ shows the simulation results for the different loss- 
cone angles of the injected electrons a^. If the injected electrons 
have the isotropic ring-like distribution (a^ - 0), the emission in- 
tensity does not exhibit a decrease at large values of resci the con- 
version efficiency of the particle energy flux into waves reaches 
12.4%. The highest emission intensity (with the conversion effi- 
ciency of up to ~ 14.5%) is reached for the horseshoe-like distri- 
butions with the loss-cone angles of about 60° - 90°. The values 
of Oc > 90° correspond to the beam-like distributions. As the in- 
jected electron beam becomes more collimated (i.e., with an in- 
creasing a^), the efficiency of the electron-cyclotron maser and 
hence the emission intensity rapidly decrease; for ffc = 120°, the 
maximum conversion efficiency of the particle energy flux into 
waves is about 7.0%. 

Finally, Fig. |5f demonstrates the effect of varying the dis- 
persion of the injected electrons in momentum Apt,. The nearly 
monoenergetic electron beams provide higher growth rate of the 
waves and possess more free energy than the beams with a larger 
dispersion in energy. Therefore, a decrease of the parameter Ap\, 
increases the efficiency of the electron-cyclotron maser One can 
see in Fig. |5]F that in the models with a smaller Ap\,, the emis- 
sion can be efficiently produced in smaller sources, with a lower 
particle escape timescale Tesc and hence by the electron beams 
with a lower density. The maximum conversion efficiency of the 
particle energy flux into waves increases with a decreasing Ap\, 



and reaches 22.4%, 14.5%, and 9.2% for A/Jb/pb = 0.1, 0.2, and 
0.3, respectively. 



4. Comparison with observations 

In this Section, we compare the results of our numerical simu- 
lations with observations. Namely, we consider the auroral kilo- 
metric radio emissions of the Earth and Saturn (for which in situ 
observations within the source regions are available), and the 
radio emission of the ultracool dwarf TVLM 513-46546 (here- 
after TVLM 513). Table [T] summarizes the main characteristics 
of the considered emission sources: the electron cyclotron fre- 
quency vb, the L-shell number of the magnetic field line where 
the source is located, the distance from the planet/star center R 
(relative to the planet or star radius Rq), the spectral intensity of 
the emission /obs at the distance d from the source; the transverse 
size R±, the electron density and energy E\, (where known), 
and the radiation energy flux per unit volume (5W/5f)obs- The 
latter parameter was estimated in the following way: since the 
emission frequency v almost coincides with the local electron 
cyclotron frequency, the observed emission in the frequency 
range dv is produced in the volume dV ^ R\dz, where the height 
range dz - (dv/v)LB and Lb is the magnetic field inhomogeneity 
scale along the magnetic field line (for the considered sources. 
Lb - R/3). Therefore, the average radiation energy flux per unit 
volume can be estimated as 



^/obs 



/obsAQ^-, 



(17) 



where AQ is the emission solid angle. We used the value of 
AQ ^ 0.22 which corresponds to the emission perpendicular to 
the magnetic field with the beam width of A^rad - 2°. 

Table |2] lists the parameters used in the simulations, as well 
as the simulation results. The effective longitudinal source size 
R. was calculated as the distance along the magnetic field line 
between the given point (i.e., where the emission with the given 
frequency is produced) and the lower (high-frequency) boundary 
of the radio-emitting region. For the auroral radio emissions of 
the Earth and Saturn, this boundary co rresponds to the electron 
cyclotron frequency of about 800 kHz (IZarkall 19981 : iLamv et aP 
2011); for TVLM 513, the high-frequency boundary of the radio 
spectrum is unknown and therefore we assumed that the radio- 
emitting region extends down to the stellar photosphere (i.e., to 
R/Rq - 1). The characteristic particle escape time from the emis- 
sion source region Tesc was calculated using Eq. and the loss- 
cone boundary was calculated using the transverse adiabatic 
invariant: 



sin ffc - vb/vb„ 



(18) 



where VBmax is the electron cyclotron frequency at the above- 
mentioned lower boundary of the radio-emitting region. The 
electron injection rate was chosen so as to obtain the specified 
particle density We in a quasi-stationary state. In all simulations, 
we assumed that the initial effective temperature of the electro- 
magnetic waves To = 10^ K. The simulation results presented 
include the radiation energy flux from unit volume (c?W/c?f)iad in 
a quasi-stationary state, the corresponding conversion efficiency 
of the particle energy flux into waves (cJW/5f)rad/('5W/50mj, and 
the emission beam width (at the 1/e level) A^rad- 

Auroral kilometric radiation of the Earth is the best studied 
among the planetary radio emissions, wit h a lot of remote and in 
situ observations (see, e.g., the reviews o f lZarkall998l:lTreumannl 
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Table 1. Parameters of the considered sources of electron-cyclotron maser radiation. 



Object 


vb 


L 


R/Ro 


lobs 


d 


R^,km 


He, cm 


£■(,, keV 




Earth (#1)° 


400 kHz 


~ 20 


- 1.61 


> 5 X 10^2' W m-- Hz^' 


1 AU 


< 100 


< 1 


- 1 - 10 


> 2.9 X 10"" 


Earth (#2)* 


~ 400 kHz 


-20 


- 1.61 


< 10-** W m-- Hz^' 


- 1100 km 


- 350 


0.2-0.5 


< 10 


< 2.5 X 10^* 


Saturn'' 


10 kHz 


-20 


- 4.75 


< 10"* Hz ' 


1 AU 


- 1000 


0.003-0.01 


< 10 


< 5.2 X 10^'^ 


TVLM513'' 


4.5 GHz 


2.15 


1.07 


< 20 mjy 


10.5 pc 


- 1000 (?) 


? 


? 


< 8.2 X 10^' 



References. IZarkal ( fl998h ; lTreumaniJ j2006h : "'' lErgun et ai] i l2000h : ''^' iLamv et all d20ia.l20lTh : '"^ iKuznetsov et al.l ( 120121) . 



Table 2. Simulation parameters and the resulting radio emission characteristics in a quasi-stationary state for the different radiation 
sources. 



Object 


vb 


R„km 


Rj_, km 


Eb, keV 






fie, cm ^ 




^md' '^'■g cm 


-3 


KJK, 




Earth (#1) 


400 kHz 


2100 


100 


10 


45° 


0.1 


0.5 


0.036 


7.9 X 10" 


-9 


0.036 


- 1.7° 


Earth (#2) 


400 kHz 


2100 


350 


10 


45° 


0.1 


0.5 


0.036 


2.4 X 10" 


-8 


0.11 


-2.0° 


Saturn 


10 kHz 


200000 


1000 


10 


5° 


0.1 


0.01 


3.4 


6.1 X 10^ 


12 


0.13 


- 1.9° 


TVLM513 


4.5 GHz 


4900 


1000 


10 


60° 


0.2 


4.2 X 10^ 


0.084 


1.0 X 10" 


-2 


0.13 


-2.8° 



0.15 




Fig. 6. Simulated quasi-stationary electron distribution in the 
magnetosphere of the Earth (model #1, R± = 100 km), a) 2D 
distribution function f(Pz,p±)- The p~ axis is reversed accord- 
ing to the conventions adopted for the northern hemisphere of 
the Earth, so that p~ < corresponds to the upward direction, 
b) ID distribution function in the energy space (integrated over 
the pitch angle). The dotted line shows the distribution of the 
injected electrons (i.e., as if the wave-particle interaction was 
absent). The simulation parameters are given in Table|2] 



0.15 




Fig. 7. Same as in Fig. |6] for the larger radiation source in the 
magnetosphere of the Earth (model #2, Rj^ = 350 km). 



l2006h . The first model in Table |2]represents a typical source re- 
gion with the transverse size of ~ 100 km, average emission 
intensity (Za rka 1998), and typical electron density and energy 
( lTreumannl2006l) . The spectrum of the terrestrial radio emission 
covers a broad frequency range; we chose vb = 400 kHz as a 



typical example. The simulated electron distribution in a quasi- 
stationary stat43 is shown in Fig. |6] One can see that the distri- 
bution is relatively weakly affected by the wave-particle interac- 
tion and does not differ significantly from the distribution of the 
injected electrons; the surface plot of the distribution function 
looks similar to that shown in Fig.[TJ). The particles are concen- 
trated in a rel atively narr ow energy range, which agrees with the 
observations (iTreu mann 2006). The conversion efficiency of the 
particle energy flux into waves amounts to only a few percent, 
but it is sufficient to provide the observed emission intensity (ac- 
tually, the source size and the electron density and energy may 
be even smaller than in the considered example). 

The second model of the source region of the terrestrial 
auroral radio emission in Table |2] corresponds to a particular 
event observed by the FAST satellite and reported by Ergu n et aH 
(2000). This event is characterized by an unusually large source 
size (R^ ~ 350 km) and very high emission intensity (far ex- 
ceeding the average values). We simulate this event using the 
same parameters as in the previous model, but with an increased 
transverse source size. As a result, the electron distribution in 
a quasi-stationary state is now more strongly relaxed (see Fig. 
|7]l and the conversion efficiency of the particle energy flux into 
waves reaches 11%. Due to an increased maser efficiency (and 
larger source volume), the emission intensity is now consider- 
ably higher than in the previous model. We have found that by 
using the observed source dimensions and electron energy and 
density, we are able to reproduce the observed emission inten- 
sity. The calculated elec tron distribution a lso looks similar to 
those observed by FAST ("Ergun et al.ll2000l) . 

Recently, the Cassini spacecraft cr ossed ffie source region 
of Saturnian kilometri c radio emission jLamv et al.ll2010i 1201 ll; 
ISchippers et al.l 1201 ih ; the basic source parameters measured 
during ffiis encounter are presented in Table[T] These parameters, 
however, seem to be rather untypical because, firstly, they cor- 
respond to the low-fre quency edge of the Saturnian radio emis- 
sion spectrum (IZarkal[l9 98). The observations were performed 
at a relatively large altitude; as a result, the loss-cone feature is 
almost absent {a^ ~ 5°) and the particle escape timescale Tesc 
is much longer than that at the Earth. Secondly, the observed 



' Note that the conversion from the 2D to the ID distribution function 
shown in Figs. l6l9l includes multiplication by p^. 
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3x10"^ 




Fig. 8. Same as in Fig.|6] for the magnetosphere of Saturn. 




Fig. 9. Same as in Fig.|6l for the magnetosphere of an ultracool 
dwarf. 



emission intensity far exceeds the average values (IZarkalll998h . 
Nevertheless, by using the model based on the observed source 
dimensions and electron energy and density (see Table|2]i, we are 
able to reproduce the observed emission intensity. The simulated 
quasi-stationary electron distribution (see Fig. [Hll agrees with 
the Cassini observations (Lamyetal. 2010; Mutel et al. 2010). 
This distribution is slightly more relaxed than that for the sec- 
ond model of the terrestrial auroral radio source; the conversion 
efficiency of the particle energy flux into waves reaches 13%. 
We remind, however, that the considered case seems to be rather 
uncommon; in a more typical source of Saturnian auroral ra- 
dio emission, the maser efficiency is expected to be much lower 
(similar to that at the Earth). 

The nearby M9 dwarf TVLM 513 represents a "classic" ex- 
ample of a radio-emitting ultracool dwarf. To date, this object 
has been dete cted as a radio sou rce in the 1.4 - 8.5 GHz fre- 
quency range dJaeger et al.ll201 ll) . although the actual spectrum 
may be broader. In addition to a quiescent component (which 
may be attributed to the gyrosynchrotron radiation), the radio 
emission includes short periodic pulses with high brightness 
temperature and ~ 100% circular polarization; the p ulse period 
coinc i des with the stella r rotation period (Hallinan et al.ll2006l 
120071; iDovle et al.1 120101: iKuznetsov etaP |2012). Most hkely, 
the periodic pulses are produced due to the electron-cyclotron 
maser instability , whic h requires a magnetic field of > 3000 G. 
IKuznetsov et al.l (1201 2h analyzed the radio light curves of TVLM 
513 and identified the likely location of the emission source 
within the stellar magnetosphere (see Table [U by assuming a 
dipole magnetic field geometry. We assume that the transverse 
size of the emission source is similar to that at Saturn (~ 1000 
km), since the radii of the ultracool dwarf and the planet are 
comparable. We also assume that the typical electron energy in 
the magnetosphere of the ultracool dwarf is similar to that at 
the Earth and Saturn (~ 10 keV), although it may actually be 
higher as suggested by the presence of a considerable gyrosyn- 
chrotron component in the radio emission. The electron injec- 



tion rate (5ne/50inj = 5 x 10^ cm"-' s"' was chosen so as to 
provide the observed emission intensity. The resulting electron 
density in a quasi-stationary state is «e - 4.2 x 10^ cm"-'. This is 
much higher than in the planetary magnetospheres, but the cor- 
responding plasma-to-cyclotron frequency ratio (wp/wB - lO""*) 
is even lower than in the sources of auroral radio emissions of 
the Earth and Saturn, and is comparable with the estimations 
for Jupiter (Lecacheux et al. 1991; Melrose & Dulk 1991). The 
simulated quasi-stationary electron distribution (see Figs. [TJ' and 
|9]i is strongly relaxed and nearly flat. Therefore we expect the 
electron distributions in the magnetospheres of ultracool dwarfs 
to differ significantly from the shell- or horseshoe-like distribu- 
tions that are typical of the source regions of planetary auro- 
ral radio emissions. Instead, we expect them to look similar to 
Maxwellian or kappa distributions, which are only slightly dis- 
torted by the parallel electric field and magnetic mirroring (i.e., 
the particle scattering on the waves is strong enough to keep the 
electron distribution close to an equilibrium state). Nevertheless, 
even these small deviations from an equilibrium distribution 
seem to be sufficient to produce an intense radio emission. 

The typical duration of radio pulses from TVLM 513 is 
about tens of seconds; moreover, the individual pulses seem to be 
caused by the star's rotation when a narrow radio b eam sweeps 
periodically over an obse rver (Hallinan et al. 2006, 20071 120081 : 
iBerger et alj|2009l : IKuznetsov et al.ii2012i) . while the parameters 
of the emission source region are stable at the timescales com- 
parable to or exceeding the rotation period (^ 2 hours). These 
timescales far exceed the estimated timescale of the particle es- 
cape from the emission source (resc ~ 0.084 s). Therefore the 
assumption that the electron distribution has reached a quasi- 
stationary state is well justified. 

5. Conclusion 

We have performed numerical simulations of the electron-cyc- 
lotron maser instability in a very low-density plasma. The used 
kinetic model included the relaxation of the electron distribution 
due to the wave-particle interaction, the injection of the energetic 
electrons with an unstable distribution into the emission source 
region, and the escape of the electrons from that region. A finite 
amplification time of the electromagnetic waves (caused by their 
escape from the source) was taken into account. The injected 
electrons were assumed to have the horseshoe-like distribution. 
We have found that: 

- The produced radio emission corresponds to the fundamen- 
tal extraordinary mode, with the frequency slightly below the 
electron cyclotron frequency and the propagation direction 
nearly perpendicular to the magnetic field. This conclusion is 
valid both for the cases when the electron distribution is sim- 
ilar to that of the injected electrons, and when it is strongly 
relaxed. The intensity of other wave modes is negligible. In 
a relatively short time (which is determined by the particle 
escape time from the emission source region), the electron 
distribution and hence the emission intensity and spectrum 
reach a quasi-stationary state. 

- Under the conditions typical of the sources of terrestrial and 
Saturnian auroral radio emissions, the dominant factor af- 
fecting the electron distribution is the particle escape from 
the emission source region. As a result, the electron distri- 
bution in a quasi-stationary state is weakly relaxed, i.e., it 
does not differ significantly from the horseshoe-like distri- 
bution of the injected electrons. The conversion efficiency of 
the particle energy flux into waves is typically a few percent. 
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although it may be higher in particular events. The emission 
escape from the source region has a "stabihzing" effect on 
the electron distribution, i.e., it reduces the maser efficiency 
and allows the particles to have highly nonequilibrium but 
long-hving (quasi-stationary) distributions. The simulated 
emission intensities and electron distributions agree well 
with the observations. 
- Radio emission from ultracool dwarfs is much more intense 
than the planetary radio emissions. Therefore we expect that 
in the magnetospheres of ultracool dwarfs, the dominant fac- 
tor affecting the electron distribution is the wave -particle in- 
teraction. As a result, the electron distribution in a quasi- 
stationary state is strongly relaxed and nearly flat; the conver- 
sion efficiency of the particle energy flux into waves is about 
10%. The radiation directivity pattern slightly differs from 
that for a shell- or horseshoe-driven instabiUty; in particular, 
the emission is directed slightly downwards. The energetic 
electrons with a relatively low density (wp/wB ~ 10"^) and 
energy of ~ 10 keV are able to provide the observed emis- 
sion intensity. 
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